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We establish the equivalence of 2D contour dynamics to 
the dispersionless limit of the integrable Toda hierarchy con- 
strained by a string equation. Remarkably, the same hierar- 
chy underlies 2D quantum gravity. 
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1. Laplacian growth. Contour dynamics takes place 
in many physical processes, where an interface moves 
between two immiscible phases. The key example of 
interface dynamics to illustrate the main result of this 
work is the Laplacian growth (LG) This process 
is dissipative, unstable, ubiquitous (applications range 
from oil/gas recovery to tumor growth), and universal: 
a steady self-similar pattern appears governed by scaling 
laws, most of which still yet to be derived |^,|^ . 

In this paper we show that an arbitrary interface dy- 
namics has an integrable structure which is the same 
as the one that underlies models of 2D quantum grav- 
ity. This structure links the interface dynamics, and es- 
pecially LG, with other branches of theoretical physics, 
where scaling laws are also expected p7| |. 

2. To be specific, we will speak about Hele-Shaw flow 
: a viscous fluid (oil) and a non- viscous fluid (water) 

are confined in a narrow gap between two parallel plates. 
The interior water domain, Dj^, is surrounded by an exte- 
rior oil domain, D- , occupying the rest of the plane. Wa- 
ter is supplied from the origin and pushes the oil/ water 
interface, C{t). Both liquids are incompressible, so oil 
is extracted at infinity at the same rate q as water is 
supplied. 

- the normal velocity of the interface is Vn = —dnP 
(the D'Arcy law); the pressure p is kept constant 
{p = 0) inside the water domain D+ (t) and on the 
interface (surface tension and viscosity of the water 
are neglected); and pressure is a harmonic function, 
V^p = 0, inside the oil domain D^{t), while p 
— (g/27r) log \Jx^ at infinity. 

This (idealized) LG problem has an important prop- 



erty: The harmonic moments of the oil domain Cfe = 
/d (t) '^~''dx dy {k — 1,2,.. ., and z = x + iy)[] do not 
change in time, while the area of water domain, grows 
linearly in time The proof: 

dC^^f Y^^I (^pOnZ-^-z-'^dnp)dC, 
at Jc{t) z Jc(t) 

because Vn — —dnP and p — along the C{t) . By virtue 
of the Gauss' theorem, it equals 

= / W {pWz^'' - z^''\7p)dxdy = -qSko- 

This property may be used as the definition of the ideal- 
ized LG problem: 

- To find the form of the domain whose area increases 
while all harmonic moments remain fixed. 

This problem is known to be ill-defined ^ . For almost 
all sets of harmonic moments, the boundary develops 
cusp-like singularities in finite time (area) . Once a sin- 
gularity occurs, the idealized LG model is no longer valid. 
Surface tension, omitted above, stabilizes the growth and 
simultaneously ruins the conservation of harmonic mo- 
ments. Simulations and experiments show that differ- 
ent mechanisms of regularization of singularities (surface 
tension, lattice etc.) exhibit the same self-similar pattern 
||l|,|| . This suggests a fixed point (or points) in the space 
of harmonic moments, which correspond to observed sta- 
ble patterns. To identify the fixed points and their scaling 
properties is the challenge of the growth phenomena. 

Approach to a fixed point requires a change of all mo- 
ments. This is the question we address in this paper. 

We present the set of differential equations which de- 
scribe the evolution of the domain under a variation of all 
harmonic moments. This prompts to a connection with 
the inverse potential problem W : to restore the shape of 



^Ci and C2 are finite: the divergence as — > 00 cancels by 
integration over argz. 
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a body from a given Newtonian potential of a uniform 
mass distribution inside the body. 

It remains to be seen whether these equations help 
to describe the pattern of growth; however, they reveal 
the integrable structure of the growth problems. We 
will show that the equations describing the evolution of 
a domain form an integrable hierarchy. Moreover, the 
very same hierarchy emerges in c = 1 string theory and 
topological gravity Q, and in 2- matrix models It 
is the dispersionless limit of the 2D Toda hierarchy |^] 
constrained by the so-called string equation ||] . 

3. To proceed further we need some known facts about 
the Schwarz function. (See e.g., jo)). An equation for a 
curve, Fc{x,y) = Fc{^^, = 0, can be resolved (at 
least locally) with respect to one of the complex vari- 
ables, say z — X — iy. The result, z = S{z), is called the 
Schwarz function of the curve C (see e.g., [^): 

(a) S{z) is a unitary operation: S{S{z)) — z; 

(b) The unit vector tangential to the curve is dz/dl — 
{dz{z)/dz)-^/^ ^ 



dzl4dld!l={dz{z)ldz)-'^l'' = l/^/Sl= VSg] 
(c) For simple analytic curves, the Schwarz function can 
be analytically continued to some strip-like domain con- 
taining the curve. 

The function S{z) can be decomposed into a sum of 
two functions S'-^^z) that are regular in D±: S{z) — 
S'(+)(z) + S'(-'(^). Under the condition S'^~\oo) = 
this decomposition is unique. The functions S^^\z) can 
be represented by a Taylor series convergent near the 
origin (which is assumed to be in -D+) and near infinity 
in D.: S(+Hz) = EZoSkz", S(-\z) = S.^z-\ 
The coefficients S±k are nothing but harmonic moments 
of the exterior, Z?_, and the interior, D^, domains: 



C±k = T 



2^'^dxdy= <t> 

Jc(t) 



S{z) dz 



2iz 



±k 



— TrS± 



fc-i 



(1) 



In other words, ^'^•'(z) is the gradient of the Newto- 
nian potential created by matter uniformly distributed 
in the interior (exterior) of C In these terms the ideal- 
ized LG problem implies that S^'^'^ does not vary in time 
and ttS-i = (Area of -D+) grows linearly in time. 

The Schwarz function is closely related to conformal 
maps. Let (j){x, y) be the function harmonically conjugate 
to 2'Kp{x,y)/q. Then w — e~^^P/«+*"^ univalently maps 
the oil domain to the exterior of the unit circle. This 
map sends w = cx) to 2 = cx). Let us write 



z{t, w) = r{t) w + Uk{t) w 



(2) 



k=0 



for the inverse map, where r is chosen to be real, so 
the map z{w) is unique. This map and the map to the 
complex conjugate domain D_, 



z{t, w 



(3) 



resolve the unitary condition for the Schwarz function 
and give it the following interpretation. If w is the image 
of a point z, then S{z) is the complex conjugate pre- 
image of w~^: S{z) = z{'w~^ (z)). 

4. The idealized LG problem has an instructive form 
in terms of the Schwarz function: 



dtS = -d^logw. 

TT 



(4) 



To derive (^) (following jl^), we differentiate z{t, w ^) = 
S(t,z{t,w)). We get dtz{t,w-^) = dtS{t,z) + 
dzS{t, z)dtz(t,w) and, by virtue of (b), Vn = Im(z(Z;) — 
St/{2iy/Sz). From log w(z) = -2ttp {x , y) / q + icf) {x , y) 
and p = const along C(t), we conclude that —dnP — 
qWz/i^iriwy/Sz). Since Vn — —dnP, we obtain (^). From 
now on we set q = tt by making a proper time rescaling. 

The equation (^) written in terms of the conformal 
maps (1), (|) has the form 0: 



{z{t,w), z{t,w ^)} = 1, 



(5) 



where we define the Poisson bracket by {/, g} = 
w{dw f dtO—dt f dyjg) . On comparing powers of w for both 
sides of (|^), we get a set of equations for the coefficients 
UkjUk of the Laurent series, (||) and (^), with fixed Ck- 
Eq. (H) suggests the Hamiltonian structure of the prob- 
lem: z{t,w), z(t,w~^) and logw, t are canonical pairs. 
5. Consider the function fl{z) defined on the curve by 



n(z) 



+ 2iA{z), zeC, 



(6) 



where we have separated the real and imaginary parts. 
Here A{z) is the area of the sector enclosed by C and 
bounded by the ray arg z and some fixed reference ray. 
Just like the Schwarz function, this function can be an- 
alytically continued within a strip containing the curve. 
Indeed, writing m) in terms of S{z), we have 



n{z) = zS{z)/2 + 2i J[S{z')dz'-z'dS{z')]/4i = J S{z')dz' 

00 00 

= J2Ckz''/TT+t\ogz-vo/2-J2c^kZ-''/7r, (7) 



fc=i 



k=l 



where vq does not depend on 2;. As is seen from the 
above definition, Q is defined up to a purely imaginary 
z-independent term. We fix it by requiring vq to be real. 

The function Q is the generating function of the canon- 
ical transformation (logui, t) (z, z) . Indeed, from (|^) 
we have S{z) = dzfl{z) and, by virtue of (^), log it; — 9fil. 
Therefore, the differential dil — S{z) dz+\ogw dt encodes 
the LG equations (^ or (||). 

6. Now we extend the differential dfl to include varia- 
tions of all higher moments. For fc > 1, let us denote 



A;=0 



^ Ck „ dfliz) 
tk = —r, tlk — —^7 — 
tt/c atk 



dn{z) 

dtk 



(8) 
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Then the multi-time Hamiltonian system is defined as 

oo 

dfl = S{z) dz + logwdt + "^^{Hkdtk — Hkdtk) ■ 

k=l 

Thus, the flows with respect to the "times" tk are 

dkz = {z, Hk} , d^z = {z, Hk} , (9) 

where dt = d/dtk, dj. = d/dtk, and {,} is the canoni- 
cal Poisson bracket introduced above. These equations 
are consistent due the symmetry relations diHk{z) — 
dkHi{z) which follow from (||). In terms of these con- 
ditions have the form of the zero-curvature equations: 

dk Hi{w) - di Hkiw) = {Hk (w), Hi (w)} . (10) 

We now proceed to calculate the Hamiltonians. Below 
we will prove that in addition to (^), for z on the curve, 
the Hamiltonians can be equivalently defined as 

Hk = -dkn{z). (11) 

(The derivative is taken at fixed z.) Then ^ gives us 

oo 

Hk = z''- dkVo/2 - dkvi z-'/l = (12) 
1=1 

oc 

= dkVo/2 + Y,dkViz~'/l, (13) 
1=1 

where we set Vk — C-fc/vr. Eq.(|lJ) implies that the Lau- 
rent expansion of the H^ at w = docs not contain 
powers of w higher than w'^. Moreover, all non- negative 
powers of w come from the first two terms of Eq. ( [l^ ) . In 
turn, Eqs. ( [l3| ) and (|^) imply that Hk does not contain 
negative powers of w. Altogether they mean that Hk is 
a polynomial in w of degree k. It reads 

Hk{w) = {z'{w))+ + ^{z''{w))o. (14) 

The symbol (/(w))+ means a truncated Laurent series, 
where only terms with positive powers of w are kept; 
(/(u)))o is the constant (z/;°) part of the series. 

It remains to prove Eq. (pT|). We first notice that 
djRei}{z) = if z belongs to the curve. This prop- 
erty is proved by differentiating the real part of (|^), 
il{z) + d{z) = |zp. The analytic continuation away from 
the curve gives il{z) + 0.{S{z)) = zS{z). Taking the par- 
tial derivative with respect to tj and restricting the result 
to the curve again, we get: 

dk^iz) + dkCiiz) + dkS{z) d-M^) = zdkS{z). (15) 

But the r.h.s. and the last term in the l.h.s. of (|l^) are 
equal since z = S{z) = dztt{z). Thus ([l5| ) reads 

dk[n{z) + n{z)] ^ 0, (16) 



where z belongs to the curve. Eq. (pT| ) follows by virtue 
of (H). Now we see that Hk and Hk (defined as in (||)) 
are indeed complex conjugates on the curve. 

Eqs. (|) or alternatively ^) together with (^, (|), 
and (^) provide an algorithm generating equations for the 
coefficients of the conformal map (||). Two first Hamil- 
tonians are Hi = rw + uo/2, and H2 — r'^w^ + 2ruQW + 
rui + Wq/2. The first equation of the hierarchy is 

dliip = dt exp{dt(p), (17) 

where = exp{dt(p)- One can see from 9fJ7 = logw that 
if is the constant term in the expansion (^: ip = vq. 

7. The unitarity condition (a) for the Schwarz func- 
tion and the properties (^ , (|l6|) of the generating func- 
tion O (which actually follow from the unitarity) impose 
important relations among the harmonic moments of any 
smooth simply connected domain and the harmonic mo- 
ments of its complement. First, from (^ one can derive 
the following sum rules: 

^ fcffcWfc =^ fctfcUfc , vi^tti+'^ktkVk-i . (18) 

k>l k>l k>2 

Second, there are symmetry relations for derivatives of 
the harmonic moments Vk = C-k/i^ of the interior do- 
main Dj^ with respect to the (rescaled) harmonic mo- 
ments tj of the exterior domain: djVk — dkVj ,djVk = 
dkVj, dtVk = dkVQ. The proof: it follows from ( |l^ ) that 
§^ HjdHk = for all j, k, then it is easy to see that 

djVk — <j> z'^dHj /2Tri — <j) z^dHk/2ni — dkVj . 

This implies that for each analytic curve C{t,tj) there 
exists a real function (prepotential) F(t, tj, tj) such that 

vj = djF , Vj ^d]F, vo^ dtF . (19) 

This function determines Hk{z) via (12). 

8. Equations (||), ( [T^ ) are familiar in the soliton lit- 
erature as the dispersionless limit |^] of the 2D Toda hi- 
erarchy. Dispersionless hierarchies of this kind Jill are 
extensions of the integrable equations of hydrodynamic 
type ||l^ to the multidimensional case. Many special so- 
lutions were found in p3[. Eq. (^ is known as the string 
equation ||^. (See also [|14[ for more recent developments 
in the 2D Toda hierarchy.) To make the contact we now 
review, following Q , the standard setup of the 2D Toda 
hierarchy and its dispersionless limit. 

The 2D Toda hierarchy is usually introduced by means 
of two difference Lax operators: 

00 

L^r{t)e''^^ +J2uk{t)e-'"'^\ 

k=l 

00 

L = r{t-n)e-^^' -f ^i2fe(i)e~'=''^', (20) 
fc=i 
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where r, Uk and Uk are functions of t and of two sets of 
independent parameters tk,tk, fc > 0. These functions 
obey the Lax-Sato equations: 



hdi^L^[Hk,L], 



where Hk = {L'') + + [L%/2, 
where Ht = {L^)- + {L%I2. (21) 



The symbol {L^)± means the part of the operator that 
consists of positive (negative) powers of the shift operator 
e''^*, and (i'^)o is the part that does not contain the 
shift operator. The first equation of the hierarchy is the 
famihar 2D Toda equation: 



(22) 



where r'^{t) — e'^^*+''^ It is also customary to con- 

sider the Orlov-Shulman operators |15| 



k=l k=2 

oo oo 

M = ^khL^ +t+^VkL-^, (23) 



fc=i 



fc=2 



where W/t — Vk{t,tj,tj), Vk = Vk{t,tj,tj) are functions 
such that the operators obey the conditions [i,^/] = 
hL, [L, M] = —TiL. These operators satisfy the follow- 
ing linear equations: L^t — z'^ , dk"^ — H^'^f, Tizdz'^ = 
Mvl/, and similarly for the bar-operators acting on 'I'. 

9. One particular solution of the 2D Toda hierarchy 
describes 2-matrix models. Consider the integral over 
two hermitian N x N matrices 



dMdM 



(24) 



(the partition function) . It has been shown that this inte- 
gral is the T-function for a special solution of the 2D Toda 
hierarchy with h = 1/N The solution is selected by 
the string equation [L, L] — h. The coefficients Vk of the 
Af -operator are given by Vk = dk logr =< Af*^ >, where 
< • • • > means an average over matrices with the weight 
( p^ ) . A scaling behavior of a proper large N limit of the 
matrix model is expected to describe 2D gravity JTzt . 

The dispersionless hierarchy is obtained in the limit 
?i — > 0. In this limit, the shift operator e'*'^' is replaced by 
a classical variable w, the Lax operators are substituted 
by their eigenvalues L — > z{w) and L —> z(w~^), and 
the operator, L~^M, becomes a function S{z). At the 
same time all commutators are replaced by the Poisson 
brackets with the symplectic structure {w,t} = w, so 
[L,L] = h turns into (I). Eq. (|l^) is the ?i limit of 
Toda's equation ( ^ ) . The wave function is replaced by 
e^/'*, where fl is the generating function of the canonical 
transformation {\ogw, t) -~f (z, z). At last, the r-function 



IS T 



0, where the function F is the 



10. To summarize, comparing the semiclassical limit 
of the Toda Eqs.( pO[), ( |2l|) , and the string equation, with 
Eqs.(||), (||), (P), (|14|), (^ of an arbitrary interface dy- 
namics, we find an exact equivalence between them. This 



is the main result of this work. The sum rules (18) for the 



harmonic moments are nothing else but (a part of) the 
PF-constraints for the r-function. It also follows from 
above that the interface dynamics is equivalent to the 
N ^ oo planar limit of the matrix model (p^): the (log- 
arithm of) the partition function of the latter is the pre- 
potential function F (p^). 

It is tempting to understand what do the robust scaling 
behavior observed in variety of growth problems 1^ and a 
scaling behavior of 2D gravity [0, (is), have in common. 
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